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ABSTTf.CT 

Conputer simulations of the differential absorption lidar technique in a space 

craft for the purpose of tenperatuie and humidity profiling indicate: 

1) Current technology ^plied to O2 and H2O lines in tte . T to .8pm wavelength 
band gives sufficiently high signal- to-noise ratios (ip to 5o for a single 
pulse pair) if backscattering b^ aerosol particles is high, i.e. profiling 
accurate to 2 K for tenperature and lo % for humidity should be feasible 
within the turbid loiter troposphere in 1 km layers and with an averaging 
over spproxiiiately ioo pirlses. 

2 ) The impact of short tenn fluctuate orjs in aerosol particle concentration is 
too big for a one laser system. Only a two laser system firing at a time lag 
of about 1 millisecond can surraount these difficulties. 

3 ) The finite width of the laser line and the quasi-random shift of this line 
introduce tolerable, partly systematic errors. 


1 . INIHODUCTION 


The retrieval of tenperature and humidity profiles from radiances measured with 
current space borne passive radianeters mainly suffers from a poor hei^t reso- 
lution. Therefore any system increasing the nurrber of independent infonrations 
- equivalent to enhanced height resolution - at the same or even slightly reduced 
accuracy level ■'s highly desirable. The present study simulates the backscatter- 
ing of laser li0:it by the atmosphere to a space craft flying at 25o or 800 km 
above the earth surface. The main goals of this study, toiperature and humidity 
profiles, need at least two wavelengths, vrtiile the aerosol backscattering profi- 
les in clear atanospheres and above clouds as well as cloud top detection can be 
determined from a single wavelength lidar. 

We report on tenperature and humidity profiles fron a two wavelength differen- 
tial absorption lidar (DIAL) . Quite a number of parameters, vrtiich all influence 
tlie backscattering of laser li^t and the accuracy of its detection have been 
set to realistic values or varied within realistic ranges in order to assess 
their inpact cn the signal-to-noise ratio (SNR) : 

- tropospheric and stratospheric aerosol particles 

- absorber amount (for humidity profiles only) 

- laser pewer (remaining within technologically feasible ranges) 

- inscability of the laser wavelength 

- width of the laser line 

- shift of the O2 or H2O line center with atmospheric pressure 

- time integration (equivalent to layer thickness desired) 

- line shsqpe of the atmospheric absorption line 

- space averaging (necessary to increase SNR) 
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2. BASIC EQUATIONS 

We start vdth the lidar equation describing the pcwer i^ceived at the detector 
F(r) after backscattering of laser power Pq by a substance at distanc r. 

P(T) s A*" exp (-2 j err'J (1) 

with Pq = emitted laser power accounting for losses within the laser optics 
yjo = correction for beam overlapping 
yjl = optical efficiency of the receiving system 
if = effective receiver area 
t = pulse length 
c = velocity of li^t 

6(r) = backscattering coefficient at distance r 
e(i) = extinction coefficient at distance r'= r 

Discretization of the integral in (1) into m finite layers with thickness Ar and 
and formulation for two wavelengths "X-i and X -2 leads to the ratio of signals 

nTr) = rMr 

fi ('T, ) JJ (- lt('r;x,)Ayj 

TJ (-2e(r'^X,i)Ar) 

u -f 

vrtiere r'= r - at/ 2 has been introduced. Writing (2) in terms of the layer i-1 
leads to a recursion formula, which - after the neglection of the wavelength 
dependence of B (certairily valid for/i\-.2nm) - is reduced to 


V('T,Xi') _ V(‘f-A'r ,X,) exp ( ~ 2. k (y'X^) Ar-) 

" P(r-4v-,X-2) exp(-2 

Resolved for the differential absorption coefficient A.k(r^ = k(r“iXi) - kCr^X^) 
we get 


A k(rr') 


- — - tn 

2AT 


. ?Cr-Ar X^)\ 

P(-r^ Xj) ■V('r - X^) j 


(3) 


If one laser wavelength, X 2 , is in a window between absorption lines, the 
differential absorption coefficient Ak(rO = k(r',Xi). New adopting the Lorentz 
line sh£^ the absorption coefficient at wavenunter v reads 


vrtiere S = line intensity depending on temperature T 
Vq = central wavenuntoer at standard conditions 
(X = line halfVddth oepending on T and pressure p 
t - line shift depending on T and p 

we introduce one of the desired parameters: temperature T. 

Accounting for the temperature and pressure dependence of S and 
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we have to handle some more parameters, 

Sq = line strength at reference tenperature 
Eq = energy of the lower state 
ICq si Boltzmann§ constant 

= molecular constant depending on degrees of freedom for i-otation 
(^ = 1 for linear molecules, = 1.5 for H2O) 

<x.Q = halfVidth at Tq,Po 
P o = 1 atm pressure 

but we also find the second desired parameter: the nuntoer density n(rO of the 
molecule having a line at Vq. Insertion of (^!) for ^=Vq and of (5) into (3) 
gives the basic equation for tarperature determination as long as n/p is known. 
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( 6 ) 


For given T and p this expression can be I'eformulated to allow for instance 
water vapor determination. 
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H,o 


Ic^ 1.132 UeC.p 
^0 Po 




(7) 


A more thorou^ treatment of the DIAL technique should account for the Doppler 
arjd Voigt line sh^ too. However, since the subsequent error equations are also 
valid for the Doppler line shape and deviations fron the Lorentz line shape in 
the troposphere are well within the errons of laser line shape knowledge, we 
omit in accordance with a statement by Megie et al (1983) a thorou^ treatment 
of Voigt line profiles. 


3. E3®0R ASSESSMENT 

We define the error of a function f(xi;i=l,. ..,m) by 


. Z. I'TTrl 
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which resolved for AXj i?eads 
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Applying ( 9 ) to ( 3 ) v;e get for the absorption coefficient 

_ i„ ( * 2 1 — 


Ak(^'; = 


ZAr r('r',X,) 


^P(rX) 


+ 2 


2 AY PCy'Xj.') 




do) 


AY . SNR, 


ay snRj. 


The approximation uses P/AP = SNR and neglects the first term in (lo), 
which is of the order of lo“^ to lo”8 for an error of lo in in the detemination 
of Ar = 1 km and for an optical depthtT’= k(r) Ar between lo ~3 and lo”l, while 
the second and the third vary from lo“l to lo "3 for SNR between 1 and loo. 
Calculating AT/T according to (9) using (6) yields 
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The application of ( 9 ) to ( 7 ) finally gives the relative eiror 
water vapor profiling. 


nH2o/^H20 for 
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4 . THE SIGNAL-TO-miSE RATIO 


Ihe error equations (lo) - ( 12 ) mainly depend on SNR at the detector. The most 
iiiportant noise sources for this study (fixed to the .75 /<ni wavelength region) 
are: 

1) shot noise caused by statistical fluctuations of photons constituting the 
sigial current 

2) background noise due to backscattered or reflected sun light 

All other noise sources have been shown to be negligible (Endemann, 1983). For 
the above noise sources rfegie and Menzies (1980) propose 
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with the number of backscattered photone received per unit time fron layer k 


N. 


2 AT 


(14) 
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and the nuntoer of background photons per unit time 
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Pj^ = signal power from layer k after passage through the front filter in front 
of the detector 
filter transmission 

■»?r = quantum efficiency of the detector 
f = frequency ( f = c/k) 

Ib = background radiance 
jT = field of view of the receiver optics 
B = filter bandwidth 
h = Planck’s constant 


5. RESULTS 

5.1 Inpact of Variations in Aerosol Particle Concentration 

Before we discuss the signal-to-noise ratio of the proposed system the errors 
due to random fluctuations of aerosol particle concentration in both the hori- 
zontal and vertical direction between two pulses of a multi-wavelength lidar 
are presented. Table 1 clearly demonstrates, that there is no way to get accurate 
profiles with only one laser shifting wavelength within o.l second. If particle 
concentration varies by only o.l percent - still unrealistically lew - between 
two pulses, the relative error calcxilated from loo trials with a randan nunfcer 
generator reaches at the maximum o.7, o.4, and 2.1 percent for temperature T, 
pressure p and number density r^j20 troposphere 

Table 1; Maximum error in percent caused by variations in aerosol particle 
concentration between two pulses. 


Variation 

AT/T 

A p/p 

An/n (H2O) 

level 

1 % 

2o.o 

8.0 

37.0 

troposphere 

o 

l.b 

2.3 

11.0 

stratosphere 


0.7 

0.4 

2.1 

troposphere 


5.2 Sigial-to-Noise Ratio (SNR) 

This section sticks to a two laser system, as shewn to be necessary within the 
previous section. The separation of pulses has to be belew a millisecond in order 
to neglect particle ccaicentration changes. The SNR calculations include the major 
noise sources: shot noise and background noise of the sunlit atmosphere as ex- 
pressvXi in (1?) - (15). Other noise sources have been handled by Endanann (1983) 
and shown to be negligible for the proposed system whose performance will be 
tested under laboratory contitions. 

SNR depends on such a variety of parameters that part of them has been fixed in 
order to assess the influence of some specifically. In accordance with the envi- 
saged hardware development we have set: '>?T=o.45, •»?r=o»l8 (Endemann, 1983a), 

A* = 0.25 m2,H= o.5 mrad, B = 1 rm, Ar = 1 km; Ib varies from 0 (nif^t) to 5 
(standard illumination) and 2oo W mr2 sr“l m“l (strong sn glint or thick clouds). 
The central wavenunber v = 12975 cnrl, Po = 25 W, and the hei^t of the space 
craft varies from 25o to 8oo km. The line parameters used are suimarized in 
table 2. The first atmosphere adopted is a standard atmosiAiere for 6o"N, July, 
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vdth 23 tan horizontal visibility at the ground and an aerosol scale height of 
1.25 km. The SNR values displayed in figures 1 and 2 for 25o and 800 tan clearly 
demonstrate : 

SNR i 10 for a single shot is only achievable vd.thin the troposphere, SNR for the 
on-line shot is drastically reduced for the stronger O2 line. However, this does 
not mean enhanced teirperature error, because the conccxnittant increase in optical 
depthtT= Ar k(r') more than conpensates this reduction. An additional constraint 
is; measurements must be taken outside sun glint areas. Within these first figures 
we did not tune aerosol backscattering in order to reach higher SNR, 


Table 2; Lines and windows chosen 


■ 

wave length X 

wavenunfcer 

— 

absorber 
and use 

+ 

line strength Sq 

halfwidth 

^0 

lower state 
energy Eq 


/Am 

cm“i 


cm 

cm-l 

cm~l 


0.7707129 

12975 • 000 

windcw; T 





0.7711989 

12966.825 

O2 ; T 

0.734 10-26 

o.c42 

1803.18 


0,7698983 

12988.728 

O2 ; T 

0.431 10-25 

o.o 42 

1422. 5o2 

M 

0.7866600 

12711.972 

H2O ; n 

0.362 10-24 

0.0655 

7o4.2l4 

h- 

0.7867201 j 

12711.000 

windcw; n 


4 



^at 296 K and I0I3.25 hPa 
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Figure 



1 : 


rower [wj 


yi^ic^ profiles of backscattered power ( ) and SNR for a standard 

afanosj^re at 60 N, July, with 23 km horizontal visibility, aerosol 
scale hei^t 1.25 km, space craft at 25o km. The thicker c’jrves stand 

rZ thinner ones for the stronger line.Pulse energy 

for all the simulations is kept at o.5 J. ^ 
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Figure 2; As fig.l, however space craft at 800 km 

Since figures 1 and 2 show SNR for one pulse pair, an averaging over loo pulse 
pairs, which will be tolerable from the point of view of horizontal resolution, 
gives SNR > loo in large parts of the troposphere for a space craft at 25o km. 

Using a standard aerosol model (Lenbble and Brogniez, 1983) for continental 
areas in the troposphere and a young ( 1 month) stratospheric aerosol layer of 
volcanic origin SNR is enhanced conpared to the atmospheric model in figures 1 
and 2, now surmounting the value io even for a space craft at 800 km under nif^it 
conditions. All the other aerosol models proposed, which were irput to our pro- 
gramne but are not shown, do not give stron^Ly different results. They all point 
to the drastic dependence on aerosol particle concentration. Rayleigh scattering 
can not give high enou^ SNR. 

5.3 Relative Temperature and Humidity Error 

The error equations (11) and (12) were esqjlored under a variety of situations. 
However, the only figure displayed alreaidy contains the main result. Die taiper- 
ature error AT/T in the lowest troposphere (figure 4) for Ar=2 km is belcw 2 K 
only for SNR2 (window wavelength) - 2oo and rather low total optical depth r“ of 
tte line. For water vapor profiling veiy similar results exist. However, the rela- 
tive error An/n at SNR2 - 2oo is only below o.l at a Eq = 9oo cnrl. Die dependence 
on t*is uncharged. 
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Figure 3* Vertical profiles of backseat tered pewer ( and SNR for a standard 

atmosphere at 6o*N, July, with continental aerosol in the troposphere 
and volcanic aerosol in the stratosphere, one month after eruption. 

Thick lines for a space craft at 25o km, thin lines for a space craft 
at 800 km. 

6. JIffiORS CAUSED BY THE FINITE WIDTH OF THE LASER LINE 

Most simulations hitherto reported neglect the laser line width. According to 
( 1983 ) present lasers have a line width of 3 pm, coxresponding to 
0.05 cm-i at 0.7/im. Korb and Weng (1982) assume o.o2 cm“l. Thus the laser line 
width is catp^ble to the halfWidth of the O 2 and H 2 O lines chosen. Since the 
above simulations assume the laser line to be monochronatic and to meet the ab- 
sorption line center we have to assess the errors thus introduced. The j^timal 
af^roach would be; revised analytic solutions. Here we adopt an inferior approach. 
We aver^ transmissions due to gaseous absorption within one layer spectrally 
synmetric to the line center. The differential absorpticai coefficients thus ob- 
tained ^ handled as before with the monochromatic ^proach. The result is a 
systematic error. Assuming a rectangular laser line of width o.o5 curl the relative 
tenperature error varies from +1.1 percent in the ipper tmosphere to -.53 percent 
w the lower troposiiiere. Assimng a triangular distribution, where the laser line 
intensity has dropped to 1/8 at o.o25 cm“i distance from the center, AT/T lies 
between -o.ol and +0.12 percent. The relative himidity error ^n/n then varies 
between +2.c and -o.95 percent. 
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SNR^, 

Figure 4 : Relative tenperature error 4T/T as a function of the signal-to-noise 
ratio for the vd.ndow position, SNR2, for different total optical 
depth as indicated. 

7. FLUCTUATIONS OF LASER FREQUENCY 

The instability of the laser frequency also causes errors, vAiich in this case will 
not be systematic. After Browell et al (1983) contenporary lasera show fluctuati- 
ons of * 0.3 F«n, approximately a tenth of the halfVddth of the absorption lines. 

We simulated these fluctuations by shifting the laser line in two different ways: 

1) random shift with equal probability for all deviations within ±0.3 pm; 

2) random shift with a Gaussian distribution and a standard deviation 6’= 0.3 pn 
and accounting for a finite triangular laser line width of 3 jan (lo? level). 

AT/T varies frcxn nearly o percent to o.l percent for both cases depending on 
hei^t level. Ihe values An/n reach 1.5 percent. This, hcwever has to be conpared 
to the An/n from (12) and thus is not intolerably hi^. All nunbers reported are 
valid for the troposphere only and after an averaging over loo realisations. 

8. INFLUENCE OF A POSSIBLE ABSORPTION LINE SHEPT WITH PRESSURE 

The abso^tion lines of atmospheric gases are not fixed to a distinct wavelength. ^ 
Their shift is negligibly small for many explications. The shift t-o^/2.75 
is cooparable to the laser line shift. Granges ftxm o.oo2 cm”l at 50 hPa and 
22o K to 0.0I8 cm~i at the ground at 29o K. A laser signal thus centered at the 
a^orption line at the tropopause is off the center to the shortwave wing of the 
line in the lower trcposphere. The systematic error reaches -0.8 percent for aT/T 
in the lower troposphere, while for the values An/n we have to expect \jp to -12 
percent. There should again be pointed out, that the r chosen is a theoretical 
estimate which needs verification. 
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9. CONCLUSIONS 

Ihe major result of this study is: 

Tropospheric taiperature and humidity profiles in 1 to 2 km layers may be deter- 
mined by a space borne DIAL using present technology to an accuracy of better 
than 2 K or lo percent respectively under the follcu'/ing conditions: 

1) Two laser system firing vdth a time lag below 1 millisecond. 

2) Averaging over at least loo pulses for a space craft at 25o km height. 

3) Averaging over some hundred pulses for a space craft at 8oo km. 

i|) Highest accuracy is found for both tenperature and humidity with O 2 and Ii20 
lines respectively at total optical depth at the line center within the 
0.1 to 0.5 range. 

5) Ihe optical depth of the pure Raylei^ scattering atmosphere in window areas 
has to be strongly enlarg^ by aerosol particles, i.e. only turbid atmospheres 
give reliable results. 

Maj’or gaps found: 

1) The inpact of laser line width and laser line position fluctuations or accur- 
acy has to be investigated in more detail. 

2) The poor knowledge of absorption line shift with atmospheric pressure hanpers 
an assessment of its inportance for DIAL frcmi space. 
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